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ABSTRACT 

In  this  report  a  numerical  study  has  been  carried  out  into  the  fatigue  life  enhancement 
of  a  through  cracked  plate  specimen  using  stop  drilled  holes  at  the  end  of  each  crack 
tip.  Also  the  effect  of  cold  working  and  the  use  of  interference  plugs  are  considered. 
While  it  has  been  found  that  under  a  static  remote  tension  load  the  use  of  interference 
plugs  gives  no  additional  static  strength,  results  from  cyclic  loads  indicate  that  the 
accumulation  of  strain  energy  per  cycle  is  less  with  interference  fit  plugs.  Also  the 
introduction  of  stationary  cracks  of  crack-length  0.066  to  1.0mm  have  been  considered 
with  constraints  to  prevent  crack  closure.  The  crack-tip  opening  displacements  of  re¬ 
initiated  cracks  at  the  edge  of  the  stop  hole  have  been  determined  using  the  finite 
element  method,  and  it  is  found  that  the  equivalent  stress  intensity  factors  are 
significantly  higher  than  predicted  by  linear  elastic  fracture  mechanics.  Furthermore 
the  use  of  these  results  together  with  the  FASTRAN  II  computer  program  predicts 
fatigue  lifetimes  which  are  in  reasonable  agreement  with  experimental  data. 
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Analysis  of  Fatigue  Crack  Growth  from  Cold- 
expanded/ interference  Fitted  Stop  Drilled  Holes 


Executive  Summary 

When  fatigue  cracks  are  found  in  aircraft  structures  the  safe  life  of  the  structure  is  of 
concern.  If  the  fatigue  life  can  be  safely  extended  then  aircraft  operating  costs  can  be 
lowered.  Previous  work  has  found  that  plastic  expansion  of  a  hole  resulting  in  residual 
compression  upon  elastic  recovery  can  result  in  a  decrease  of  the  mean  stress  under 
cyclic  loads  resulting  both  in  an  increase  of  life  and  an  increased  critical  crack  length.  In 
an  experimental  program  recently  carried  out  at  AMRL  a  procedure  was  used  to  stop 
drill  the  crack,  cold  work  the  stop  hole  and  use  interference  fit  plugs.  This  resulted  in  an 
increase  of  stress  level  by  a  factor  of  2  for  fatigue  life  of  10,000  Falstaff  blocks. 

In  this  report  a  theoretical  investigation  is  carried  out  into  the  fatigue  life  enhancement  of 
a  cracked  plate  specimen  using  a  combination  of  stop  drilling,  cold  expansion  and 
interference  fitting.  Firstly  the  stress/ strain  state  in  the  area  in  which  cracks  re-initiate 
from  the  cold  worked  hole  is  analysed  using  the  ABAQUS  finite  element  program.  This 
involves  an  elastic-plastic  and  non-linear  contact  analysis.  Validation  of  this  work  has 
been  carried  out  using  closed  form  solutions  for  plugs  in  holes  for  an  infinite  plate, 
corresponding  to  (1)  cold  expansion,  (2)  removal  of  the  mandrel,  (3)  insertion  of  an 
interference  plug,  and  (4)  application  of  a  remote  load.  It  has  been  found  that  in  the 
absence  of  an  applied  remote  load  the  residual  stress  distribution  in  the  area  in  which 
cracks  would  develop  is  approximately  the  same  as  would  be  expected  in  a  cold 
expanded  hole  in  an  infinite  plate.  Also  it  has  been  found  that  under  static  tension  loads 
use  of  cold  working  and  interference  plugs  gives  no  additional  strength,  however  results 
obtained  from  cycling  of  loads  indicate  that  accumulation  of  strain  energy  per  cycle  is 
considerably  less  for  cold  worked  and  interference  fit  plugs,  and  hence  a  considerable 
improvement  in  life  would  be  expected.  Secondly,  the  introduction  of  short  stationary 
cracks  of  crack-length  of  0.066  to  1.0mm  has  been  considered  with  constraints  to  prevent 
crack  closure.  The  crack-tip  opening  displacements  of  re-initiated  cracks  at  the  edge  of 
tire  stop  hole  have  been  determined  using  the  finite  element  method,  and  it  is  found  the 
equivalent  stress  intensity  factors  are  significantly  higher  than  those  obtained  by  linear 
elastic  fracture  mechanics.  From  these  results  together  with  the  FASTRAN  II  computer 
program,  fatigue  lifetimes  are  predicted  which  are  in  reasonable  agreement  with 
experimental  data  for  these  specimens. 

Overall  it  has  been  found  that  the  combination  of  interference  fit  plugs  with  cold 
expansion  of  the  stop  drilled  holes  significantly  reduce  the  crack  driving  force  and  hence 
extends  the  fatigue  life  of  a  cracked  structure.  Results  achieved  in  this  report  indicate  that 
this  is  a  promising  procedure  for  the  life  extension  of  RAAF  aircraft  structures. 
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Symbols 


a 

a 

e 

8t 

X 

e 

c 


<*z 


p 

V 

A K 
A K„ 


eff 


material  constant 

also  denotes  the  constraint  factor 

strain 

yield  strain 

crack  tip  opening  displacement 

interference 

angle 

stress 

tangential  stress 
radial  stress 

stress  perpendicular  to  plate 
yield  stress 

radius  to  re-yielded  zone 
poisson's  ratio 
poisson's  ratio  for  mandrel 
effective  stress  intensity  factor 
threshold  stress  intensity  factor 


ai  re-initiated  crack-length 

2a0  length  of  long  crack 

C  material  constant 

CTOD  crack  tip  opening  displacement 

E  Young's  modulus  of  aluminium 

Em  Young's  modulus  of  mandrel 

dn  displacement  variable 

G(ai,r)  stress  intensity  of  a  crack  ai  due  to  point  force  at  distance 

In  variable 

J  J  integral 

K  stress  intensity  factor 

Kp  plastic  stress  intensity  factor 

c  radius  to  elastic-plastic  boundary 

m  work  hardening  exponent 

n  work  hardening  exponent 

p  interference  pressure 

R  hole  radius 

r  radius 

ux  x  displacement  near  crack  tip 

uy  y  displacement  near  crack  tip 


1.  Introduction 


As  a  result  of  high  acquisition  costs  of  aircraft,  life  enhancement  techniques  are  being 
increasingly  used  to  reduce  overall  operating  costs.  Generally  the  growth  of  fatigue 
cracks  determine  the  safe  life  of  the  structure.  Life  enhancement  techniques  involving 
removal  of  cracks  in  fastener  holes  and  fitting  interference  inserts  are  well  known. 
However  experimental  work  carried  out  by  (Finney  et  al  ,1996)  has  found  that  the 
fatigue  lifetime  of  a  cracked  plate  specimen  in  which  the  crack  tips  have  been  stop 
drilled  and  cold  expanded  is  also  considerably  enhanced.  The  geometry  of  the  cracked 
plate  is  shown  in  Figure  1.  The  crack  has  a  length  of  2 .a0  =  20 mm  and  is  referred  to  as 
the  long  crack.  Numerical  and  experimental  work  has  also  been  carried  out  by  (Vulic 
1997)  for  cracks  emanating  from  a  single  hole  in  which  a  mandrel  is  used  to  cold  work 
the  hole.  However  the  aim  of  this  report  is  to  understand  the  mechanism  by  which  the 
fatigue  life  is  enhanced.  An  investigation  is  carried  out  into  the  stress  state  when  the 
hole  is:  (1)  cold  expanded,  (2)  the  mandrel  is  removed,  (3)  the  insertion  of  an 
interference  plug  and  (4)  subject  to  a  remote  applied  load.  The  interference  (A)  is 
defined  as  the  ratio  of  the  difference  in  hole  diameters  between  the  hole  and  mandrel 
to  the  diameter  of  the  hole.  In  the  work  carried  out  here  the  interference  level  for  cold 
expansion  is  4%  and  interference  for  the  plugs  is  .67%. 
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Figure  1.  Stop-drilled  crack  with  an  interference  insert 
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In  this  analysis  short  cracks  of  0.066  to  1.0mm  will  be  considered  in  the  area  in  which 
fatigue  cracks  are  expected  to  re-initiate.  It  is  well  known  that  when  the  area  of 
plasticity  associated  with  a  growing  crack  is  large  in  comparison  to  the  cracklength, 
linear  elastic  fracture  mechanics  (LEFM)  is  no  longer  valid.  Hence  a  more  accurate 
analysis  of  short  cracks  can  be  achieved  through  the  use  of  an  elastic-plastic  analysis. 

2.  Analysis 

The  F.E.  analysis  has  been  carried  out  using  ABAQUS/Standard.  This  program  has  the 
required  contact  analysis,  non-linear  geometry  and  plastic  stress  analysis  capability. 
The  material  is  a  2024-T851  alloy  and  is  idealised  as  a  perfectly  elastic  and  perfectly 
plastic  stress  state,  with  material  properties  as  shown  in  Table  I.  Currently  for 
plasticity,  no  closed  form  solutions  exist  for  the  interference  bolt  in  a  hole  containing  a 
crack.  However  closed  form  solutions  do  exist  for  the  case  of  the  expansion  of  a 
mandrel  in  the  hole  of  an  infinite  plate,  which  does  not  contain  a  crack. 


Table  I  Material  properties  for  2025-T851. 


Young's  modulus  (MPa) 

Poisson's  ratio 

Yield  stress  (MPa) 

Thickness 

(mm) 

72188. 

.33 

435. 

3.6 

The  following  closed  form  solutions  will  be  used  for  the  validation  of  the  analysis. 


Figure  2.  Definition  of  parameters  for  analytical  expression. 
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2.1  Step  1:  cold-expansion  of  stop-hole 

For  the  case  of  hole  in  an  infinite  plate  (see  Figure  2),  the  stress  distribution  in  the  plate 
and  the  plug  can  be  determined  analytically  (Timoshenko  and  Goodier,  1970). 
Assuming  that  the  plate  material  obeys  an  elastic-perfectly  plastic  stress-strain 
relationship,  the  stresses  in  the  plate  after  cold-expansion  by  a  mandrel  are  given  by 
(Jost,  1988), 


r  if  r  ~\  crs  if  r  ]  a,  1  r 

—  =  —j=  21n— -1  ,  —  =  -;=  21n-  +  l  ,  —  =  -7=21n- 

r0  V3L  c  J  cr0  V3  L  c  J  <r0  V3  c 


in  the  plastic  region  surrounding  the  hole  and 

t,  '"is1 (f) '  ■ <r- =°  <2> 

in  the  elastic  region.  Here  c  denotes  the  boundary  between  plastic  region  and  elastic 
region,  which  can  be  determined  from  the  following  equation. 


where  E,  v,  <?o  are  respectively  the  Young's  modulus,  Poisson's  ratio,  yield  stress  of  the 
plate  material,  and  Em  and  vm  are  the  Young's  modulus  and  Poisson's  ratio  of  the 
mandrel. 

2.2  Step  2:  removal  of  expansion  mandrel 

Upon  the  removal  of  the  mandrel,  recovery  can  either  be  elastic  or  elastic-plastic 
depending  whether  or  not  re-yielding  occurs.  The  condition  for  re-yielding  is 

-~>4e&  1.65  (4) 

K 


When  the  recovery  is  entirely  elastic,  the  stress  distributions  are 


cr  i  f  r  "If  c  JR) 

a-o  V3|L  c  J  [  R  JVr7 

cre  1  ff  r  1  f  c 

a0  S  L  c  J  L  R  j \r) 


—  =  -t=<  21n-  +  l 


o-o  V3[L  c  j  L  R  . 

for  R  <r  <c  and 

cr  1  \(  cV  f  c  Y  R 
—  =  — 7=^  -  -  21n— +  1  — 


f-ifel'-hri'l'l 


for  r  >  c . 
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When  re-yielding  occurs,  the  radius  (p)  of  the  re-yield  region  is  given  by 


and  the  stresses  in  the  re-yielded  zone  are: 


For  the  zone  R<r  <  p: 
a r  _  _  2  R  <7g_  _ 

&0  S  r'aQ 
For  the  zone  p<r  <c: 


For  die  zone  r  >  c : 


(10) 

(11) 

(12) 


2.3  Step  3:  insertion  of  interference  plug 

For  small  interference,  the  deformation  around  the  hole  should  remain  elastic, 
especially  when  the  material  around  the  hole  is  already  under  compression  as  a  result 
of  cold-expansion. 

The  interference  pressure  is  given  by, 

EX 

P  -  £  (14) 

(1  +  v/)  +  (1  +  vJ(1-2vJ  — 

m 

and  the  corresponding  elastic  stress  distribution  is  given  by 


The  final  stress  distribution  can  be  readily  obtained  by  superimposing  the  solutions 
determined  for  the  above  three  steps.  Comparison  with  finite  element  results  will  be 
discussed  in  the  next  section  together  with  the  results  corresponding  to  cold-expanded 
stop-hole  fitted  with  an  interference  plug. 
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3.  Finite  Element  Analysis  of  Cold-expansion  of  a 

Stop-hole 

3.1  Step  1:  cold-expansion 

The  stress  distributions  around  a  stop-hole  determined  by  an  finite  element  analysis 
are  shown  in  Figure  3  together  with  the  stress  distribution  corresponding  to  a  hole  in 
an  infinite  plate.  It  can  be  seen  that  in  the  region  6  <  100°,  there  is  a  reasonable 
agreement  (with  a  difference  about  10%)  between  the  stress  distributions  (both 
tangential  and  radial  stresses)  in  an  infinite  plate  and  a  stop-hole,  indicating  that  the 
influence  of  the  long  crack  is  mainly  limited  to  a  region  near  the  intersection  between 
the  crack  and  the  hole. 


Angular  coordinate  0 


Figure  3.  Stress  distributions  around  a  cold-expanded  hole  and  stop-hole  at  4%  expansion. 

3.2  Step  2:  removal  of  mandrel 

During  the  removal  of  the  mandrel,  further  compression  yielding  occurred  in  the 
tangential  direction  around  the  hole,  as  marked  by  the  lower  compressive  tangential 
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stress  shown  in  Figure  4  as  compared  to  that  shown  in  Figure  3.  Again,  in  the  region 
0  <  100° ,  the  stress  distribution  is  nearly  unaffected  by  the  presence  of  the  long  crack. 


Angular  coordinate  0 
Figure  4.  Residual  tangential  stress  after  the  removal  of  mandrel. 

3.3  Step  3:  Insertion  of  interference  plug 

Upon  insertion  of  an  interference  plug  into  the  cold-expanded  stop-hole,  the  residual 
tangential  stress  around  the  hole  is  reduced  in  magnitude;  the  finite  element  results  are 
shown  in  Figure  5  together  with  the  results  for  the  analogous  problem  of  a  hole  in  an 
infinite  plate.  It  is  seen  that  except  near  9  =  0  for  the  two  cases,  there  is  a  considerable 
difference  between  the  tangential  stress  distribution  around  the  hole.  Consequently  the 
analytical  method  developed  for  a  hole  in  an  infinite  plate  is  no  longer  applicable  and 
finite  element  analysis  has  to  be  adopted.  The  tangential  stress  distribution  away  from 
the  hole  edge  ( 6  =  0 )  is  shown  in  Figure  6,  indicating  a  good  correlation  between  the 
theoretical  prediction  and  the  finite  element  results  in  the  case  of  a  hole  in  an  infinite 
plate.  It  is  also  important  to  note  that  as  far  as  the  tangential  stress  distribution  away 
from  the  hole  is  concerned,  the  analytical  predictions  are  still  in  good  agreement  with 
the  finite  element  results. 
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Angular  coordinate  0 


Figure  5.  Stress  distribution  after  the  insertion  of  an  interference  plug  (X-  0.67%  ). 


Distance  from  hole  edge  r/R 

Figure  6.  Distribution  of  the  tangential  stress  around  of  a  cold-expanded  hole  in  an  infinite 
plate  followed  by  insertion  of  a  interference  plug  (X  =  0.67% ). 
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Figure  7.  Tangential  stress  at  edge  of  cold  expanded  hole  versus  remote  stress  for  case  of  no 
plug  and  interference  plug. 

4.  Application  of  Remote  Load  to  Combined 
Interference  Fitting  and  Cold-Expansion  of  Stop-holes 

To  examine  the  overall  benefit  of  combined  interference  fitting  and  cold-expansion  of 
the  stop-hole  in  improving  fatigue  life  of  cracked  structures,  remote  tension  is  applied 
to  the  complete  geometry.  The  tangential  stresses  ahead  of  the  hole  crack  at  different 
applied  load  levels  are  shown  in  Figure  7.  For  comparison  purposes,  the  tangential 
stress  distribution  corresponding  to  the  case  of  cold-expanded  stop-hole  without  an 
interference  plug  is  also  included  in  the  figure.  The  Young's  modulus  assumed  for  the 
plug  is  2 lOOOOMPa,  and  poisson's  ratio  is  .3  .  One  important  observation  that  can  be 
made  from  the  results  is  shown  in  Figure  7.  The  rate  of  increase  in  the  tangential  stress 
at  the  hole  edge  for  the  case  fitted  with  an  interference  plug  is  only  slightly  lower  than 
the  case  without  interference  plug.  This  implies  that  under  tension  dominated  loading, 
the  interference  plug  may  not  provide  any  extra  benefit  than  a  cold-expanded  stop- 
hole.  A  comparison  of  the  tangential  stress  at  0=0  at  the  edge  of  the  hole  versus  the 
remote  applied  load  shown  in  Figure  8,  shows  no  improvement  for  the  interference  fit 
plug.  However,  it  is  expected  that  when  the  far  field  load  has  compression 
components,  as  normally  encountered  in  many  aircraft  structures,  the  addition  of  an 
interference  plug  is  beneficial  in  eliminating  local  reverse  yielding  near  the  stop-hole. 
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This  is  further  illustrated  in  Figure  9  which  shows  the  deformed  geometry 
corresponding  to  a  remote  stress  of  200  MPa.  Clearly,  due  to  the  presence  of  the  long 
crack,  the  stop-hole  opens  up  and  the  plug  is  totally  detached  from  the  hole.  The 
contact  modelling  procedure  involves  a  master-slave  concept  in  which  the  mandrel 
corresponds  to  the  master  and  the  aluminium  plate  the  slave.  The  solution  is  iterative 
in  that  the  initial  contact  surface  must  be  assumed  and  an  improvement  on  this 
assumption  is  obtained  after  each  solution,  until  convergence  is  achieved.  The  change 
in  geometry  may  require  a  large  displacement  analysis  to  be  carried  out.  In  the  next 
section  the  advantages  of  interference  will  be  considered  under  remote  cycling  loads. 


0  1  2  3  4  5 

Distance  away  from  hole  edge  r/R 

Figure  8.  Tangential  stress  distribution  ahead  a  cold-expanded  hole  fitted  with  an  interference 
plug  subjected  to  remote  load. 
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t 


Long  crack 


Re-initiated  crack 


Figure  9.  F.E.  mesh ,  Vi  symmetry,  deformed  plot  corresponding  to  remote  stress  of  200MPa, 
shows  that  the  interference  plug  has  no  influence  for  maximum  tensile  loads. 

5.  Stress  Analysis:  Interference  Plug  Subject  to  Cyclic 

Load 

5.1  Tangential  stress  and  load  level 

In  order  to  understand  the  fatigue  behaviour  of  a  short  crack  an  analysis  has  been 
carried  out  for  the  first  5  cycles  where  the  applied  stress  ranges  from  0-200MPa.  For 
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this  configuration  only  a  long  crack  exists,  with  prior  expansion  of  the  hole.  The  stress 
to  be  considered  is  the  tangential  stress  located  at  the  point  at  which  a  short  crack 
would  seem  most  likely  to  initiate.  Two  cases  are  considered,  the  first  corresponds  to 
the  case  in  which  a  plug  is  not  used  during  loading,  while  the  second  case  involves  the 
use  of  an  interference  plug.  In  Figure  10  a  plot  of  tangential  stress  versus  applied  stress 
is  shown  for  the  case  of  no  plug.  For  this  case  the  stress  cycle  starts  at  an  applied 
remote  stress  of  zero,  as  shown  by  the  arrow.  The  tangential  stress  at  this  point  is  - 
500MPa  due  to  cold  expansion.  Note  that  the  stress  at  which  yielding  occurs  for  plane 

strain  is  cr0 .2/^3  =  500 MPa  .  Increasing  the  applied  load  stress  to  lOOMPa  results  in 
a  linear  variation  of  tangential  stress  to  a  maximum  value  of  500MPa  at  which  yielding 
occurs.  Increasing  the  applied  stress  to  200MPa  maintains  the  tangential  stress  at 
500MPa.  Unloading,  as  shown  by  the  arrow,  to  an  applied  stress  of  80MPa  results  in  a 
linear  variation  of  tangential  stress  to  approximately  -500MPa.  Unloading  to  zero 
maintains  the  tangential  stress  at  approximately  -500.MPa.  Subsequent  cycles  follow 
this  identical  path. 


Figure  10.  Variation  of  tangential  stress  with  applied  stress  for  case  in  which  no  plug  is  used 
during  loading. 


At  this  point  it  is  useful  to  make  a  comparison  with  the  case  in  which  an  interference 
plug  is  present  during  loading.  The  results  for  this  case  are  shown  in  Figure  11,  where 
the  arrows  show  the  load  cycle.  In  comparison  to  Figure  10  significant  differences 
exist.  Firstly,  the  presence  of  an  interference  plug  results  in  an  initial  stress  of  -200MPa 
at  zero  load,  in  comparison  with  -500MPa  shown  in  Figure  12.  In  Figure  11  the  initial 
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loading  to  160MPa  shows  a  non-linear  response,  which  is  due  to  the  influence  of  the 
plug.  Further  loading  to  200MPa  results  in  a  tangential  stress  of  500MPa.  The 
unloading  curve  is  similar  to  that  in  Figure  10,  except  that  a  maximum  compressive 
stress  of  -550MPa  occurs.  The  reason  for  this  is  that  under  compression  the  plug  does 
not  allow  the  hole  to  deform  and  relieve  the  radial  contact  pressure  at  that  point, 
which  due  to  poissons  ratio  increases  the  tangential  stress.  Note  that  as  shown  by  the 
arrows,  the  second  cycle  does  not  follow  the  first  cycle,  however  all  subsequent  cycles 
are  identical.  In  comparison  to  the  no  plug  case,  the  most  significant  differences  occur 
with  the  gradients  of  the  compressive  loading/ unloading  curves. 


Figure  11.  Variation  of  tangential  stress  with  applied  stress  for  case  of  interference  plug  used 
during  loading. 

5.2  Stress/strain  curves 

The  location  of  interest  is  still  at  the  region  at  which  otherwise  a  small  crack  would 
initiate.  For  this  location  the  stress/ strain  plot  corresponding  to  the  no  plug  case  is 
shown  in  Figure  12.  Again,  the  arrows  show  the  initial  loading  and  cyclic  behaviour. 
These  curves  exhibit  the  perfectively  elastic  perfectively  plastic  behaviour  during 
loading  and  unloading.  Clearly  for  tangential  stresses  greater  than  500MPa  yielding  is 
occurring  and  strains  change  significantly.  Note  that  all  cycles  show  an  identical  path. 
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Figure  12.  Stress/strain  plot  for  case  of  no  plug  during  loading. 


For  the  case  of  an  interference  plug  used  during  loading,  the  results  are  shown  in 
Figure  13,  where  the  arrows  show  the  initial  loading  corresponding  to  the  lowest  value 
of  strain,  and  the  subsequent  load  path.  The  lowest  value  for  strain  occurs  at  a 
compressive  stress  of  -200MPa  which  corresponds  to  the  fitting  of  an  interference  plug. 
For  the  first  cycle  after  the  yield  stress  has  been  reached  a  considerable  strain  occurs 
during  loading.  However  all  subsequent  cycles  result  in  a  considerably  reduced  range 
of  strain  values  in  comparison  to  the  first  cycle.  The  hysteretic  plastic  work  per  cycle 
can  be  regarded  as  a  correlating  parameter  for  fatigue  damage  accumulation.  The 
comparison  of  stress/strain  curves  for  the  no  plug  case  and  interference  plug  are 
shown  in  Figure  14.  Except  for  the  first  cycle,  it  is  evident  that  the  hysteretic  plastic 
work  for  the  no  plug  case  is  approximately  double  the  case  in  which  the  interference  fit 
plug  is  used.  As  a  result  the  expected  fatigue  life  for  the  case  in  which  a  plug  is  used 
would  be  expected  to  be  significantly  greater.  Note  also  that  with  the  use  of  an 
interference  plug  a  slight  increase  in  stress  amplitude  occurs. 
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Figure  13.  Stress/strain  plot  for  case  of  use  of  interference  plug  during  cyclic  loading. 


plug  duing  loading 
no  plug  during  loading 


Figure  14.  Stress/strain  comparison  plots  for  case  of  plug  during  loading  and  case  of  no  plug 
during  loading. 
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6.  Computation  of  Stress  Intensity  Factor  of  Re¬ 
initiated  Crack 

6.1  Analytical  approach 

When  the  enhanced  specimens  are  subjected  to  remote  loading,  depending  upon  the 
level  of  the  applied  load,  plastic  deformation  may  occur  ahead  of  the  stop-hole, 
leading  to  new  crack  growth,  as  depicted  in  Figure  15.  Here  we  will  first  examine  the 
growth  behaviour  of  a  re-initiated  crack;  the  issue  of  crack  initiation  will  be  addressed 
later. 

The  prospective  stress  distributions  along  a  newly  initiated  crack  at  various  applied 
stress  levels  are  shown  in  Figure  8.  For  a  small  crack  of  length  a,  re-initiated  at  the  tip 
of  the  stop-hole,  the  stress  intensity  factor  can  be  computed  using  a  Green  function 
(Grandt,  1975;  Ball,1990) 

K  =  £'  G(r, ax )aeg{r)dr  (16) 


Figure  15.  Coordinate  and  notations  for  crack  re-initiation  from  stop-hole 
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where  the  function  G(rfl\)  denotes  the  stress  intensity  factor  of  a  crack  a\  caused  by  a 
point  force  acting  at  position  r.  The  exact  Green  function  does  not  appear  to  exist  for  a 
short  crack  emanating  from  a  stop-hole.  In  the  following  we  will  make  use  of  the 
limiting  solutions  for  ax  « R  and  ax  »  R  to  provide  lower  and  upper  bound 
solutions. 


For  ax  «R,  the  problem  can  be  treated  as  an  edge  crack  in  an  infinite  plate 
(Tada/1974)/ 


G,(r,at) 


2  lJ-OJjr/a,)3'2 
yjl  —  (r/  a,)2 


(17) 


and  for  a,  »  R,  the  Green  function  is  that  for  a  point  force  acting  on  a  line  crack  in  an 
infinite  plate,  which  is  well  known. 


Gl(r,al)  =  - 


\n{a  +  a,) 


2  a0  +ax+r 
ax-r 


Here  the  subscripts  s  and  l  are  used  to  denote  the  Green  functions  pertaining  to  the 
short  crack  bound  and  the  long  crack  bound.  As  seen  in  Figures  16  and  17  ,  these  two 
solutions  provide  reasonable  bounds  to  the  stress  intensity  factor. 


Length  of  Re-initiated  Crack  (mm) 


Figure  16.  Stress  intensity  factor  versus  crack  length  for  remote  stress  equal  to  200  MPa. 
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Length  of  re-ininitated  crack  (mm) 

Figure  17.  Stress  intensity  versus  crack  length  for  short  and  long  crack  bounds. 


6.2  Method  for  computation  of  stress  intensity  factors 

The  following  is  a  brief  account  of  the  basis  on  which  the  plastic  stress  intensities 
are  computed,  for  a  monotonic  loading.  A  more  detailed  account  is  given  in  the 
Appendix.  The  method  outlined  here  is  not  applicable  to  cycling  loading.  For  fully 
plastic  solutions,  (Goldman  and  Hutchinson,  1975)  give  the  displacement  field  at 
the  crack  tip  for  mode  I  opening  as: 

u,=a.eaK /''•"%($)  ...(19) 

where  a  is  a  material  constant 

£0  is  a  reference  value  of  strain  (yield) 
n  is  the  work  hardening  exponent 

where  Ke  is  the  plastic  strain  intensity  factor 

Finite  Element  solutions  have  been  considered  (Shih,  1981)  for  the  case  of  small 
scale  yielding.  An  evaluation  of  the  crack  tip  opening  displacement  (CTOD)  is 
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defined  at  a  point  in  which  a  45°  line  intersects  the  crack  profile  shown  in  Figure 
18.  The  equation  is: 


r  —  u 


X 


.(20) 


Figure  18.  Intersection  of  45°  line  and  the  crack  profile. 


The  value  of  S  that  satisfies  equ(20)  is  given  by: 
8=dn  —  ~  CTOD 

ao 

Since  J  =  K2  /  E  we  have: 


K  = 


Ecr0CTOD 


where  cr0  is  the  yield  stress. 


...(21) 


...(22) 


6.3  Discussion  on  analytical  and  F.E.  analysis 

It  has  been  found  that  to  compute  the  plastic  stress  intensity  factor  using  the  F.E. 
method  it  is  necessary  to  have  a  fine  mesh  with  crack  tip  elements  of  the  size  of 
.0001mm.  In  Figure  19  the  crack-face  displacements  are  shown  at  the  crack  tip  for  a 
crack-length  of  0.097mm.  The  intersecting  45°  line,  as  discussed  in  section  6.1  is  also 
shown  and  gives  the  CTOD  solution  from  which  stress  intensities  may  be  computed. 
Each  data  point  behind  the  crack  tip  represents  a  node  in  the  mesh,  and  it  is  evident 
that  in  this  case  values  of  CTOD  are  only  achieved  for  load  levels  of  lOOMPa  and 
above. 
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Figure  19.  y  displacement  versus  distance  behind  crack  tip,  for  a  range  of  applied  stresses 
corresponding  to  a  crack-length  of 0.0979mm. 

As  mentioned  previously  a  number  of  short  cracks  of  crack-lengths  of  a  =  0.066  to 
1.0mm  have  been  considered  for  a  remote  applied  stress  of  200MPa  and  are  shown  in 
Figure  20.  The  corresponding  stress  intensities  are  those  shown  in  Figure  16.  Results 
from  cycling  the  applied  remote  stress  are  shown  in  Figure  21,  where  the  arrow 
corresponding  to  a  zero  CTOD  defines  the  initial  loading  and  other  arrows  show  the 
cyclic  loading  path.  In  this  case  only  the  CTOD  has  been  considered  since  the 
procedure  for  computing  stress  intensity  factors  corresponding  to  unloading  is  more 
complex  than  the  procedure  for  monotonic  loading  outlined  in  section  6.1  .  Note  that 
except  for  the  first  cycle  of  loading  the  cycling  follows  an  identical  path. 
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Figure  20.  Values  of  CTOD  versus  cracklength  for  a  range  of  cracklengths,  enhanced, 
A.  =.67%.,  remote  load  =  200MPa 


Figure  21.  Variation  of  crack  opening  displacement  with  applied  stress  for  5  cycles, 
a,  =  0.0997 mm,  Stress  ratio=0,  a0  -  20 mm . 


It  is  clear  from  Figure  17  that  for  the  same  applied  stress  of  100  MPa,  the  driving  force 
for  a  re-initiated  crack  increases  initially  and  then  decreases  below  zero,  indicating  a 
crack  arrest  even  when  a  new  crack  has  been  initiated  at  the  hole  edge.  This  means  that 
the  enhanced  geometry  should  yield  infinite  fatigue  life  for  an  applied  stress  equal  to 
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100  MPa.  This  is  primarily  caused  by  the  compressive  tangential  stress  between  0.2 
and  0.8  times  the  hole  radius. 

When  the  applied  stress  is  increased  to  150  MPa,  as  shown  in  Figure  8,  the  tangential 
stress  becomes  positive  everywhere,  and  hence  any  new  crack  initiated  at  the  hole 
edge  will  continue  to  grow  to  failure. 

When  the  applied  stress  is  further  increased  to  200  MPa,  some  degree  of  plastic 
yielding  occurs  at  the  edge  of  the  hole.  However,  when  the  above  results  are  compared 
with  those  determined  from  finite  element,  as  indicated  by  symbols  in  Figure  16,  the 
analytical  method  considerably  underestimates  the  stress  intensity  factor.  This  is  most 
likely  due  to  the  fact  that  the  influence  of  plasticity  near  the  hole  was  not  accounted  for 
in  the  analytical  method.  The  extent  of  the  plastic  deformation  was  approximately  20% 
of  the  stop-hole  radius.  The  results  shown  in  Figure  16  suggest  that  once  the  newly 
initiated  crack  reaches  about  0.3mm,  it  will  grow  at  much  the  same  rate  as  the  un¬ 
enhanced  crack. 


6.4  Fatigue  life  prediction 

This  is  further  demonstrated  in  Figure  22,  where  the  experimental  data  are  taken  from 
(Finney  et  al,  1996).  It  is  seen  from  Figure  22  that  FASTRAN  II  is  able  to  provide  a  very 
good  prediction  of  the  fatigue  crack  growth  curve  for  the  un-enhanced  specimen. 

In  this  case  the  growth  rate  is  determined  from  the  following  equation  (used  in 
FASTRAN  II): 


da 

dN 


(23) 


where  the  constants  corresponding  to  2024-T851  may  be  determined  from  the  Damage 
Tolerant  Handbook  and  are  given  by: 

C  =  7.4x1  O'10 
A Kth  -  1.0  MPayfin 
m  =  2.93 

a  =  1.9  (constraint  factor  required  by  FASTRAN  II) 

In  the  case  of  the  enhanced  specimen,  once  the  crack  length  exceeds  0.3mm,  the 
prediction  is  in  good  agreement  with  the  experimental  data.  At  the  moment,  it  is  not 
clear  how  we  can  predict  the  crack  growth  behaviour  for  crack  lengths  less  than 
0.3mm,  which  is  the  key  to  determining  the  total  fatigue  life  of  the  enhanced  structure. 
For  short  crack  lengths,  local  plasticity  near  the  hole  edge  plays  a  significant  role  in 
affecting  fatigue  crack  growth  rate.  Work  is  currently  in  progress  to  develop  a 
methodology  which  is  capable  of  dealing  with  short  cracks  growing  from  a  notch 
plastic  field. 
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Figure  22.  Peak  tensile  stress  versus  number  of  loading  blocks  for  a  plain  cracked  plate  and  a 
cracked  plate  enhanced  by  a  combination  of  interference  fitting  and  cold-expansion  of  stop-hole. 
Symbols:  experimental  data;  solid  curves:  predictions  using  FASTRAN II. 
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7.  Conclusions 


The  following  conclusions  can  be  drawn  from  the  present  study  of  a  particular 

improved  life  enhancement  method:  combined  interference  fitting  and  cold-expansion 

of  stop-hole. 

1.  In  the  absence  of  an  applied  remote  load,  the  residual  stress  distribution  along  the 
prospective  path  of  a  newly  initiated  crack  near  a  cold-expanded  stop-hole  (with  a 
long  crack)  is  approximately  the  same  as  that  which  would  exist  around  a  cold- 
expanded  hole  in  an  infinite  plate. 

2.  Upon  the  application  of  a  remote  load,  the  interference  of  an  inserted  plug  would 
not  only  reduce  the  stress  concentration  ahead  of  the  stop-hole  but  also  reduce  the 
level  of  hoop  stress  at  the  onset  of  plastic  yielding  due  to  the  beneficial  effect  of  the 
interference  pressure  on  plastic  yielding. 

3.  For  the  case  of  cyclic  loading  it  has  been  found  that  the  combination  of  interference 
fitting  and  cold-expansion  of  a  stop-hole  would  significantly  reduce  the  magnitude 
of  hysteretic  work  around  the  stop  hole  although  some  increase  in  stress  amplitude 
is  noticed. 

4.  After  the  crack  is  re-initiated  ahead  of  the  stop  hole,  the  value  of  the  crack  driving 
force  is  found  to  be  considerably  greater  than  the  Green's  function  approach 
(LEFM). 

5.  Based  on  the  numerical  results  of  the  crack  driving  force,  which  accounts  for  notch 
plasticity,  fatigue  life  predictions  are  made  using  a  crack  growth  code.  Provided  a 
suitable  initial  crack-length  is  selected,  good  correlation  can  be  obtained  with 
experimental  values. 
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Appendix:  Method  for  Computation  of  Stress 

Intensity  Factors 


The  following  analysis  is  the  basis  on  which  the  plastic  stress  intensities  are 
computed,  as  shown  in  Figure  16,  for  a  monotonic  loading.  They  are  not,  however 
applicable  to  cycling  loading. 


Fully  plastic  solutions  are  considered,  (Goldman  and  Hutchinson,  1975),  and  the 
relationship  between  stress  and  strain  is  given  in  the  following  form: 


where 


a  is  a  material  constant 
s0  is  a  reference  value  of  strain  (yield) 
cr0  is  a  reference  value  of  stress  (yield) 
n  is  the  work  hardening  exponent 


•(1) 


For  fully  plastic  solutions,  (Goldman  and  Hutchinson,  1975)  give  the  displacement 
field  at  the  crack  tip  for  mode  I  opening  as: 


uy  =  a.s0K£rV(n+1)uy(6>) 

•~<2) 

where  Ke  is  the  plastic  strain  intensity  factor 

This  solution  is  the  same  as  that  proposed  by  (Hutchinson,  1968)  and  (Rice  and 
Rosengren,  1968)  (HRR)  for  the  case  of  small  scale  yielding.  The  HRR  near  tip 
behaviour  is  expressed  as: 

^ 5v(0,„) 

...(3) 

(*.») 

...(4) 

In  order  to  relate  crack  opening  displacement  to  /  (Tracey,  1976)  has  expressed  the 
HHR  singularity  as: 


EJ 


\  i/(«+d 


\a.a0Inr  ) 


°iMn) 


...(5) 
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Sy  = 


a.a0 


EJ 


Ka.<y0Inr) 


b/(b+ 1) 


sJe,n) 


...(6) 


where  E  is  Young's  modulus 
<r0  is  the  yield  stress 
/  is  the  /  Integral 

Finite  Element  solutions  have  been  considered  (Shih,  1981)  for  the  case  of  small 
scale  yielding.  An  evaluation  of  the  crack  tip  opening  displacement  (CTOD)  is 
defined  at  a  point  in  which  a  45°  line  intersects  the  crack  profile  shown  in  Figure 
Al.  The  equation  is: 


r  —  u 


X 


...(7) 


The  choice  of  a  45°  angle  does  have  the  advantage  of  restricting  the  CTOD  to  near 
tip  values.  The  terms  that  satisfy  equs(7)  are  given  (Tracey,  1976)  as: 


5,  /  2  =  ^- 


EJ 


N  b/(b+1) 


a.crj, 


o  •*»' 


„B+1  "y 


«v(«) 


a.an 


u  = - 

*  E 


EJ 


\a.<r0In. 


b/(b+1) 


,B+1  "y 


«v(») 


...(8) 

...(9) 


Figure  Al.  Intersection  of  45° line  and  the  crack  profile. 
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The  value  of  §  that  satisfies  equs(7),  (8)  and  (9)  is  given  by: 


St 


=  CTOD 


Since  J  =  K2  /  E  we  have: 


...(10) 


K  = 


\Ect0CTOD 

V  ~ 


...(11) 


we  also  have  from  (Tracey,  1976): 


for  small  values  of  ux  we  have: 


dn 


\  E  J 


...(12) 


...(13) 


Values  of  dn  from  (Tracey,  1976)  have  been  plotted  in  Figure  A2  corresponding  to 
plane  strain.  Hence  from  equ(ll)  the  mode  I  stress  intensity  can  be  computed  in  a 
plastic  strain  field. 
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